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AMENDED SPECIFICATION 
1 

ACOUSTIC ASSESSMENT OF THE HEART 

5 

TECHNICAL FIELD OF THE INVENTION 

The present invention relates to non-invasive cardiac and respiratory 
monitors. More specifically, diagnostic tools and methods for evaluating the 
heart and respiratory system and improving the function of the failing heart are 
10 disclosed. 

BACKGROUND OF THE INVENTION 

Cardiopulmonary monitoring of clinically unstable patients is essential 
15 for early detection of potentially life threatening changes in the functionality of 
the heart or respiration. Cardiopulmonary monitoring Is presently done on 
hospitalized patients in critical care, such as intensive care units (ICU), coronary 
care units (CCU) and peri-operatively as well as in other specialized sections of 
the hospital. Certain types of cardiopulmonary monitoring are also performed on 
20 patients outside the hospital, such as patients who suffer from asthma, high 
blood pressure, or cardiac antiythmias. Cardiopulmonary monitoring is also 
used during exercise performance tests of patidnts and athletes. Cleariy, 
monitoring and eariy detection of Impending catastrophes are highly desirable in 
clinical medicine. 



Continuous monitoring of cardiopulmonary wellbeing is presently focused on 
physiological parameters that characterize, individually, the activity of the heart 
and respiration. These parameters include monitoring of the electrocardiogram, 
blood pressure in the arterial, venous and pulmonary circulations and cardiac 

5 output for assessing the cardiovascular system. The airway pressure, 
respiratory rate, tidal volume, flow rate, pulse oximetry, exhaled CO2, and 
esophageal pressures are the respiratory parameters that are commonly 
monitored. More recently continuous overnight monitoring of breath sounds was 
introduced to detect and quantify wheezing activity in asthmatics. Additional 

10 methods were introduced to monitor patients inflicted with congestive heart 
failure (CHF) to identify early decompensation. These methods are based on 
continuous monitoring of the electrical impedance of the thorax, but have not 
generated reliable levels of sensitivity and specificity. Monitoring of respiratory 
crackles has also been used as an early sign of lung congestion. 

15 The activities of the heart and the lungs are well known to be closely 
inten'elated. The heart rate, blood pressure and blood flow into and out of the 
heart are influenced by the breathing cycle. The pressure inside the chest 
becomes more negative during inspiration to enable inflow of air into the lung 
alveoli. This sub-atmospheric pressure also affects the heart, blood vessels and 

20 blood flow. In particular, the increased intrathoracic negative pressure expands 
the right atrium and ventricle, dilates and elongates the blood vessels, amplifies 
the ventricular filling and changes the position of the inter ventricular septum. In 
addition, the changes in lung volume modify the afferent neuronal activity in the 
vagus nerve, leading to modulation of the heart rate during respiration. During 

25 quiet breathing the intrapleural pressure decreases from about -3 mm Hg to 
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about -6 mm Hg. This causes dilation of tlie intrathoracic segment of the vena 
cava, increased venous return to the right atrium and ventricle. The increased 
diastolic filling of the right ventricle amplifies its stroke volume into the pulmonary 
circulation by the well-known Starling mechanism. At the same time, the 

5 inspiratory displacement of the abdominal content by the contracting diaphragm 
increases the intra abdominal pressure, which further pushes blood into the 
thoracic cavity. The intrathoracic pressure changes during breathing also 
influence the output of the left ventricle, but to a lesser degree. During deep 
breathing, or in pulmonary diseases that affect the mechanical properties of the 

10 ainA^ays and the lung parenchyma these, respiratory swings of flow, resistance 
and volume are greatly exaggerated (Mountcastle V. O. Medical Physiology, 
12^ ed. Mosby Company, St. Louis 1968). 

Positive pressure ventilation squeezes the pulmonary capillaries and 
increases the resistance to blood flow through the lungs. These pressures may 

15 diminish the output of the right ventricle due to increased afterload, while at the 
same time reducing the output of the left ventricle due to fall in its diastolic filling 
and preload. These phenomena are known to result in wide fluctuations of the 
cardiac output and blood pressure. Similar fluctuations are also seen during 
resisted breathing, cough and isometric muscle straining. 

20 Monitoring of heart sounds is well known to medicine and physiology 

for many years, even before the invention of the stethoscope by Laennec in 
1819. The first and second heart sounds are associated with closing of the atrio- 
ventricular and ventricular outlet valves, respectively. They are loud and distinct 
sounds with somewhat different amplitude and template in different areas of 

25 precordial auscultation. The first and second heart sounds are modified, and 



sometimes are even completely missing during diseases of the heart or the 
lungs. Additionally, the 3^ and 4^^ heart sounds are well known to be 
associated with defects of the left ventricular filling during diastole. Rumbling 
sounds heard in between the first and second sounds (systole) or the second 

5 and first sounds (diastole) are appropriately called systolic and diastolic 
murmurs, respectively and are associated with abnormal blood flow through the 
narrowed or malfunctioning heart valves. Information on heart sounds, 
phonocardiography and the art of interpretation of heart sounds is readily 
available in many text books and articles, such as in "Rapid interpretation of 

10 heart sounds and murmurs" by Emanuel Stein, Abner J. Delman (Editors), 
Williams & Wilkins, 1997, 4th edition, the contents of which are incorporated 
herein by reference. 

The respiratory changes of cardiac activity are well known. In 
particular, the phenomenon of "Pulsus Paradoxus" is a recognized sign of 

15 severe asthma and ainA^ay narrowing. It is defined as a decline of greater than 
12 mm Hg (in some texts 20 mm Hg) in the systolic blood pressure during 
inspiration. Detection of Pulsus Paradoxus in a dyspneic patient is an ominous 
sign that calls for aggressive and immediate intervention. On the other hand, 
cardiac arrhythmia, the acceleration and deceleration of the heart rate during the 

20 respiratory cycle is often a benign and normal phenomenon, particularly in 
young children. The only available information on the effect of respiration on the 
heart sounds per-se is on the width of splitting of the second heart sound and on 
certain cardiac murmurs. Otherwise, no information is available on the changes 
in the heart sounds induced by respiration and on the extent of these changes 

25 during various cardiopulmonary conditions. 
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The failing heart condition is a cause for severe morbidity and mortality 
rates. This condition can be helped by a broad range of assist methods. These 
include pharmaceutical agents (e.g., digitalis and other positive inotropic 
agents), mechanical support (e.g., intra aortic counter pulsation, and 

5 intravascular coronary artery stent), and various types of electrical stimulation 
(e.g., ventricular resynchronization therapy, Guidant's Contak CD, Medtronic's 
multi electrode epicardial pacing. Impulse Dynamics' timed refractory period 
current). All heartbet synchronized methods must be monitored as regards 
their performance. This feedback can determine if, and to what extent the 

10 assistance is effective. Such monitoring is either performed in intervals (e.g., by 
periodic detenmination of the left ventricular ejection function), or continuous 
(e.g., by monitoring the left ventricular pressure - Remon Medical Technologies 
LTD of 7 Halamish Street, Caesaria Industrial Park, 38900, Israel). To be 
effective, the mechanical and electrical assist methods should be synchronized 

15 with the cardiac cycle so that they can be activated at specific phases of the 
cardiac cycle. For example, in US patent 6,285,906 the contents of which are 
incorporated herein by reference, impulse dynamics provides current stimuli at 
certain precise phases of the heart cycle. To do so, most methods use the EKG 
(electrocardiographic) signal to find the appropriate timing within the cardiac 

20 cycle for activation of the assistance. However, the EKG may not provide 
sufficiently detailed or accurate information about the timing cardiac cycle, or 
may be contaminated by electromagnetic noise. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 . is a schematic description of the signal collection and processing device 
providing synchronized features of the cardiac activity; 

5 Figure 2, is a schematic description of the entire process in accordance with a 
preferred embodiment of the invention, performing propbabilistic analysis. 

Figure 3. is a schematic description of the decision loop involving heartbet 
synchronized control in accordance with the invention. 
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DETAILED DESCRIPTION OF THE INVENTION 

The present invention provides a novel method for monitoring the interrelated 
functionality of the heart and the respiratory system by tracking the changes in 
5 cardiac sounds with respect to the activity of the heart during the respiratory 
cycle. A novelty of the new approach is that the timing of each heart contraction 
within the breathing cycle is used as a cardinal parameter in the analysis of the 
heart sounds relative to the electrocardiographic signal. By doing so, the heart 
sounds generated during similar segments of the respiratory cycle, for example, 

10 the first part of inspiration, are analysed together to generate respiratory 
synchronized stable templates and features of the heart sounds and their 
interaction with the EKG. Thus, the extent of changes in EKG - modulated heart 
sounds parameters during each phase of the respiratory cycle may be tracked 
over time to detect variability of the breathing related swings that indicate 

15 worsening of the status of the cardiopulmonary system. In addition, the invention 
facilitates tracking of the stable synchronized cardiac sounds features over time 
that are indicative by themselves of the well - being of the mechanics of the 
heart and lungs. It is further disclosed that the detemiination of the stable 
synchronized cardiac sounds features is classified by the timing and template of 

20 the EKG. Together, these new methods provide valuable information on the 
status and mechanics of the complete cardiopulmonary system. 

The invention overcomes two major hindrances to extraction of stable 
features from heart sounds: the variability of the cardiac sounds during the 
respiratory cycle and changes in the sounds caused by variability of the 

25 electrical activity of the heart. The new method extracts features such as 
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amplitude, duration, frequency content, template and chirp components of the 
heart sounds and clusters these features with respect to their timing in the 
respiratory cycle and /or their underlying EKG morphology. This method 
preserves the attributes of the cardiac acoustic signal that are othenA^ise 

5 eliminated or greatly diminished by grouping the sounds or their features without 
respect to their respiratory timing and EKG, as is currently learned by the 
existing art. The extraction of synchronized stable features greatly diminishes 
their stochastic variability, so that any changes in the deterministic components 
of the heart sounds, such as those caused by changes in the mechanical activity 

10 of the heart, are more readily and accurately detected and characterized. The 
acoustic properties of the heart sounds reflect the mechanics of the cardiac 
contraction, the integrity of the cardiac structure and dynamics of flows, volumes 
and pressures in the cardiovascular system. Therefore, detection of changes in 
the synchronised stable sounds relative to a baseline may be used to alert the 

15 patient, his or her caregivers, or an automated algorithm or device that can be 
prompted to provide a remedy to the patient's condition. This new method, while 
preserving complete non-invasiveness, provides far more comprehensive 
monitoring of the patient than the cunrent existing methods such as the EKG. 

The grouping of synchronised stable features of the heart sounds with 

20 or without respect to the EKG template defines a baseline extent of Within 
Breath Variability (WBV). Any changes of the extent of variability of the sounds 
or their features within the respiratory cycle relative to a previous period, such as 
the baseline, can be readily and accurately detected. These alterations of the 
WBV may be caused by increased pressure swings within the chest due to 

25 altered breathing pattern or changes in the resistance, compliance or other 
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mechanical properties of the lung. Such changes are often seen in asthma 
patients, patients inflicted by congestive heart failure (CHF), patients who are 
artificially ventilated, anaesthetised, or otherwise critically ill. 



5 Analysis of the heart sounds 

In Fig. 1 to which reference is now made the main stages in collecting 
and processing of the signals, in a preferred embodiment of the invention are 
shown. Cardiac acoustic signals 20, respiratory acoustic signal 22 and EKG 
electric signal 24 are collected by the appropriate detectors. The cardiac 

10 acoustic signal and the respiratory acoustic signal are separated and 
conditioned, by module 28, which includes in a preferred embodiment of the 
invention filtering, denoising, amplification, multi-sensor analysis for removal of 
background noise and recording artifacts. Feature extraction module 30 of the 
invention performs segmentation. Heart sounds are segmented into S1, S2, S3, 

15 S4, the breathing activity is segmented into the Heart sounds are segmented 
into S1, S2, S3, S4, the breathing activity is segmented into five respiratory 
stages; first half of inspiration, second half of inspiration, first half of active 
expiration, second half of active expiration and the post expiratory pause if 
present, and the EKG signal is segmented into its components; P, QRS, T 

20 waves. The segmented input data undergo timing analysis by module 32. 
Temporal segmentation feature parameters are extracted from the segmented 
data by module 34. The parameters are estimated individually from heart 
sounds appearing at the different respiratory stages. The EKG signal is also 
analyzed with respect to the respiratory stages. 
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Heart sounds are acquired via a single microphone or a plurality of 
microphones located on the chest, back and/or neck of the patient. After signal 
conditioning that includes pre-amplification, the signal is band-pass filtered so 
that the respiratory component may be separated from the heart sounds. The 

5 isolated sound structures are then detected in the heart sounds signal and 
classified as S1, the first heart sound, and S2 the second heart sound. The 
degree of isolation (which is related to the beating rate) is estimated earlier from 
the data using a density plot of the inter-beat duration. A temporal window, 
related to the time interval variability, is used to distinguish between S2 and S1. 

10 Further validation of the S1/S2 segmentation may be provided by the EKG 
signal where the timing of S1 is known to con-elate with the QRS sequence and 
the timing of S2 corresponds to the T wave. After the different sounds are 
segmented, an attribute corresponding to the respective respiratory cycle is 
attached to each of them, so that they can be further collected into groups with 

15 their corresponding timing within the respiratory cycle. Averaging the different 
sounds and/or their attributes in the corresponding respiratory cycle forms 
adaptive templates. Averaging is performed within each timing interval on the 
raw signal, the frequency content representation, and on a time-frequency 
representation obtained by short time Fourier transfomn, Wigner distribution 

20 analysis or a more adapted time-frequency representation using continuous 
wavelet transform or best basis and discrete wavelet representation. The latter 
methods include internal denoising (compression) which is achieved by 
removing the low energy coefRcients from the transformed signal representation 
before performing the inverse transfomi. 
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The time-frequency representation can be used for parameter 
estimation of tfie chirp properties (mainly S2). These properties include the 
chirp slope, namely the frequency modulation of the signal, start and stop 
frequency as well as slope, and the amplitude changes in the signal. Properties 

5 related to the time-frequency representation are extracted using conventional 
curve fitting techniques. Features related to the amplitude changes, e.g. 
increase in negative amplitude of S2, frequency content as is recorded by 
Fourier transform, and instantaneous frequency changes which is recorded in 
the time-frequency representation. The parameters corresponding to each of the 

10 above categories, and their variability may be empirically estimated from the 
data and compared to a data base or baseline values of these parameters. The 
variability together with the estimation of the mean, may further be used for the 
estimation of statistically significant changes in each of the parameters, both 
during a respiratory cycle and when compared with a normal baseline of the 

15 specific patient or of previously recorded and labeled data. 

The separation of the deterministic and stochastic signal components 
is done during entropy-based denoising, e.g. with wavelet or best basis 
representation. The wavelet coefficients with smallest variance, which 
correspond to smallest contribution to the energy of the collection of heart 

20 sounds, are set to zero. This parameter shrinkage removes the part of the signal 
that is independent of the heart sounds and can be considered noise. 

A non exhaustive list of temporal segmentation parameters and 
features, with respect to the respiratory stages extracted from the heart sounds, 
respiratory cycle and EKG useful in the cardiac assessment of the invention is 

25 provided next: 
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I . Amplitude and energy of S1 and S2 in each segment, (peak 
and RMS) where S1 and S2 are the first and second heart 
sounds, respectively. 

2. Timing of S1 and S2 relative to the EKG and relative to each 
other whereas the EKG represents the electrical activity of the 
heart, measured with skin-surface electrodes. 

3. Delta time between occurrence of S2 and S1 and occurrence of 
S1 and the S2 of the previous heart cycle. 

4. Delta time between the occurrence of S1 and occurrence of the 
QRS wave. 

5. Delta time between the occurrence of S2 and occurrence of the 
P wave. 

6. Delta time beween the occurrence of the P wave and the R 
wave. 

7. Frequency content and instantaneous frequency content of S1 
and S2. 

8. Cord length of S1 and S2. 

9. Polarity of first and largest component of S1 and S2. 

10. Amplitude incline and damping (decaying) rate constant of S1 
and S2. 

I I . Chirp parameters of SI and S2. 



12. Duration of whole cardiac cycle (reciprocal of the heart rate) as 
well as of each heart sound segment (i.e. of S1, S2, S3, S4, or 
murmurs). 

13. Signature (or shape) specific for each sensor by using signal 
processing methods and/or polynomial fit. Differentiate between 
deterministic elements (e.g. SI -4) and stochastic elements such 
as murmurs and/ or artefacts. Determine separately for each 
frequency range. Obtain template pattern of the heart sounds 
signature for each segment of the respiratory cycle. 

14. Cross correlation between various precordial sites. An 
expansion of this is creating a dynamic acoustic map of the 
anterior chest using many sensors, each displayed on an array 
with brightness indicating amplitude and color(s) indicating the 
frequency content. The dynamics of the changing map can be 
tracked for changes during the respiratory cycle. 

15. Calculate the deconvolution of the heart sounds/ ballistographic 
signal in order to estimate the mechanical activity at the source 
that generates the signs. 

16. An AM/FM Modulation decomposition of the heart sound to 
analyse independently amplitude changes and frequency 
changes in the heart sounds. 

17. Amplitude multiplied by the width of each heart sound which is 
an approximation to the temporal length of the heart sounds 
times the peak amplitude of each heart sound. Energy of some 
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components of the heart sound as is measured in a 
time/frequency representation which can be adapted to the 
signals. 

Using the above parameters and features, deterministic and stochastic 
5 components of the features and parameters are extracted. This defines the 
acceptable boundaries for each feature and parameter. Such boundaries are 
determined by the extent to which the parameters vary whilst the patient is in 
stable condition. In a preferred embodiment, these boundaries are defined 
by the mean and n x SD of each parameter and feature for each the 
10 respiratory stages, where n is 2 for example. In particular, the (potentially non- 
linear) boundaries of the hyper-parameter space of the full set of the 
parameters are determined. Thus, for example, it is possible that a patient can 
be at a point where each of the parameters and features is within the defined 
boundaries, but the full set of parameters is outside of the boundaries. 

15 A probabilistic model, which is estimated in a preferred embodiment, 

via calculation of maximum likelihood, is applied to the processed data to further 
refine the model estimation. This binds together the temporal progress from one 
respiratory stage to another (via a hidden Markov model). At this stage, 
additional information from other sources such as blood pressure, oxygen 

20 saturation, patient weight, chest electrical impedance and other means, is fused. 
The full temporal model and set of segmented features are then analyzed to 
detenmine scoring of deviation from nonnality, or from the acceptable patients' 
boundaries. 
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The probabilistic analysis assumes a Gaussian distribution to the 
continuous features of the model and a binomial (sometimes multinomial) 
distribution to the categorical parameters. Following the estimation of the mean 
and standard deviation, significance or abnormality of an observed features may 
5 be determined and using a probability model which takes into account the 
probability of a collection of features (taken from an already collected and soon 
to be collected data), one can estimate the probability of various abnormal 
events. 

In Fig. 2 to which reference is now made, the whole sequence of steps 
10 from the signal transduction to model display is described schematically. In step 
50 signals are transduced to produce a conreponding electric cun^ent. In step 52 
the signals are conditioned. In step 54 the signals are segmented to express the 
segments of the physiological rhytmicity. In step 56 Parameters are extracted, 
and In step 58 the deterministic components are separated from the stochastic 
15 components. In step 60 the probabilistic analysis is performed, accepting 
physiological data at step 62, and in step 64 scoring and display are provided. 
Although not shown, the probabilistic analysis performed reflects on the signal 
segmentation to influence the segmentation process. 

20 Sensors for collecting signals of the heart and respiratory system 

In general, apart from the ubiquitous EKG which senses electrical 
signals of the body, all other sensors taking part in the signal collection in 
accordance with the present invention are acoustic or mechano - acoustic 
sensors. Such sensors transduce mechanical pressure, or acoustic energy to 

25 electrical energy further to be analyzed electronically. Acoustic sensors for 
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sensing activity of the heart to determine cardiac mechanical performance, 
which include typically rhythmicity. fine cardiac temporal cycle variability and 
muscle contractility, can be patched on the patient's chest or can be implanted 
subcutaneously. Conditioning respiratory activity sounds and respiratory activity 
5 mainly for the determination of respiratory rhythm is preferably performed using 
impedance plethysmography. Several sensors can be used concomitantly for 
obtaining data about the sounds. 

Electronic hardware for processing the data and providing 
10 diagnostic aids and control signals 

Transduced electric signals from the sensors, as well as EKG signals, 
are collected by data collecting hardware, sampled and digitized on a dedicated 
processor or in a computer. All processing and statistical calculations and model 
assessment are done in the computerized environment. A GUI or a printer is 

15 typically used to produce a report to the medical team in charge. In that case 
that control signals are required, this can be supplied by a controller linked to 
the processor or computer wherein the processing is performed. If more than 
one sensor is used for obtaining sounds data, a procedure for correlation 
between the sensors may be implemented. The implementation of such a 

20 procedure is typically carried out for reducing noise. 

Medical uses of the invention 

A simple example of the method of the invention includes segmenting 
the breathing cycle into exhaling and inhaling segments, and recording 
25 contracting heart sound at these inhaling and the exhaling segments. Sound 
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amplitude samples are collected and then averaged, average A for the set of 
samples representing the exhaling segment period and average B for the set of 
samples representing the inhaling segment period. Statistical difference is 
calculated between the averages A and B, providing evidence regarding 
5 quantitative difference in magnitude between the two sets. The physiological 
significance of the difference is such that the larger difference is indicative of a 
pronounced tendency of the heart to be limited by the expanding lungs. 

These parameters and changes may be used to detect early changes 
associated with pulmonary and/or hemodynamic alterations. These include CHF 

10 (changes in the heart's mechanical activity, the lungs specific gravity (density) 
and increased fluctuations of intrathoracic pressure); mechanical ventilation 
(changes in PEEP, airway resistance), changes in hydration and blood volume 
(shock) (cardiac mechanical activity, increased breathing activity), hypertension 
(increased cardiac mechanical activity), myocardial infarction (changes in 

IS cardiac contraction sequence). 

The detection of each of the above-listed conditions is done by 
comparing the respiratory synchronized cardiac sounds or features thereof to a 
patient-specific baseline, or to a template consisting of data from an aggregate 
of normal and abnormal patients, A system of the invention or any system using 
20 the method of the invention are either therapeutic or diagnostic. 

Cardiac synchronization methodology and usage in diagnostic 
and therapeutic devices 
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Assessment of the mechanical condition of the heart performed in 
accordance with the invention can be used to help treating the failing heart. An 
acoustic heartbeat synchronized method comprises, in accordance with the 
present invention, one or more sensors of heart sounds that can pick up and 

5 transduce the heart sounds of a patient into electrical signal and means for 
amplifying the signal, typically but not necessarily a sensor for breathing activity 
and means for synchronizing the signals. The heart sounds timing signal can be 
used to trigger the activity of a heartbeat synchronized device, either directly or 
with a delay. Such delay may be pre-determined or modified automatically, 

10 based on information on the heart sounds amplitude, the interval between them, 
or their amplitude and frequency content. 

In Fig. 3, to which reference is now made, an automated control 
mechanism a heartbeat synchronized device using a synchronization algorithm 
of the present invention, is shown. Acoustic cardiac signals are transduced 80 

15 and respiratory signals 82 are transduced to be collected by the monitor 84. 
Subsequently, they are processed at step 86 by the method of the invention to 
extract accurate parameters of the heart. If the heart is regarded failing at 
diagnostic step 88, the heartbeat synchronized device is activated. 

The heart sounds timing signal is used to actuate a demand-type or 

20 synchronized cardiac pacemaker. This pacemaker is chronically implanted, or 
temporarily used to defibrillate a fibrillating heart, such as used in automatic or 
manually triggered electric defibrillators. Several therapeutic and diagnostic 
systems can benefit from using the heartbeat synchronization concept 
embodied in the invention. Such systems are intra-aortic balloon pump, left 

25 ventricular cardiac assist device, coronary angiography diagnostic devices. 
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cardiac imaging devices such as, but not limited to, CT, MRI and SPECT 
diagnostic technologies. 

5 Monitoring variability changes 

The grouping of synchronized stable features of the heart sounds with 
or without respect to the EKG template defines a baseline extent of within breath 
variability (WBV). Any changes of the extent of variability of the sounds or their 
features within the respiratory cycle relative to a previous period, such as the 

10 baseline, can be readily and accurately detected. These alterations of the WBV 
may be caused by increased pressure swings within the chest due to altered 
breathing pattern or changes in the resistance, compliance or other mechanical 
properties of the lung. Such changes are often seen in asthma patients, patients 
inflicted by congestive heart failure (CHF), patients who are artificially ventilated, 

15 anesthetized, or othenA/ise critically ill. When such changes are found they can 
be used to alert the patient, his or her caregivers, or an automated algorithm or 
device that can be prompted to provide a remedy to the patient's condition. 

Short duration analysis 

20 Known in the art are changes in the duration of each heartbeat in 

association with the breathing cycle. These changes are usually detected by 
measuring the distance between consecutive "R" waves in the EKG that are 
called "the R-R interval". The changes in the duration of the R-R interval over 
time are referred to as heart rRate variability" (HRV). HRV is considered as 

25 indicative of the activity of the sympathetic and parasympathetic control of the 
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heart. It is substantially reduced in gravely ill patients. Recently, changes 
(reduction) in HRV have been associated with particulate matter (PM) air 
pollution and are believed to be an important mechanism in the increased 
mortality due to PM pollution. Other conditions, such as the sleep apnea 

5 syndrome, chronic cardiac ailments and the use of certain medications are also 
known to be associated with changes in HRV. 

In accordance with the present invention, analysis of the short-term 
duration changes of systolic time (Tsys) and diastolic time (Tdias) is performed 
and studied with the respiratory cycle as well as the other amplitude and signal 

10 morphology features that have been discussed above. Of particular interest are 
correlation between short temi changes which are indicative of a well-heart 
being such as increase in Tsys which leads to a an increase in the amplitude 
and sharpness of S2, indicating an increase in duration that is due to more blood 
flow into the heart, or a negative correlation between Tsys and Tdias, indicating 

15 more elaborated breathing (possible air-way obstruction) but healthy heart 
dynamics, vs. changes in Tsys that are not correlated with Tdias and vice versa. 

The method of the invention can be used to differentiate between 
inhalation and exhalation phases of the respiratory cycle. It can also be used to 
determine the template and variations of S3, S4 and cardiac murmurs. 

20 

Use of the invention in irregular breathing instances 

The method of the invention can be applied in the course of irregular 
breathing. Forced breathing interruptions are examples in which the method of 
the invention may be used to analyze the condition of the heart. More particular 
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examples for such irregularities are breath hold after deep respiration and 
Valsalva maneuver while agonal respiration is a pathological instance of 
interrupted respiration. 

5 Use of the method of the invention in diagnosing extra cardiac 

blood vessels 

The method of the invention can be applied to improve the diagnosis 
of extra - cardiac blood vessels using MR! (magnetic resonance imaging) as 

10 well as other imaging methods. The application with MR! is known in the art as 
MRA (magnetic resonance angiography), and includes applying a magnetic field 
to the region of interest, and mapping the effect of the magnetic field on various 
constituents of the region of interest. The technique is very useful, but the 
lengthy data acquisition time of the procedure causes blurring of the image 

15 obtained. In order to improve the image, synchronizing the MR imaging with the 
heartbeat and other imaging methods or with the heartbeat as synchronized 
versus respiratory activity improves the clarity of the image. 



